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ABSTRACT 


This  paper  describes  a  finite  element  method  for  analysing  the  behaviour  of flaws  in 
thin  fibre  composite  sheets  which  are  required  with  a  bonded  overlay.  The  method  is  an 
extension  of  the  authors'  previous  work  on  the  repair  of  metal  structures.  As  illustrative 
examples  the  repair  of  a  hole  and  a  crack  in  various  graphite  I  epoxy  laminates  is  discussed 
and  it  is  shown  that  in  each  case  a  bonded  overlay  of  either  boron  I  epoxy  or  titanium  pro¬ 
vides  effective  repair. 


NOTATION 


Cartesian  coordinate  system 
Transverse  shear  stresses 
Transverse  shear  strains 
Adhesive  shear  stresses 

Thickness  of  the  sheet,  adhesive  and  patch  respectively 
Transverse  Shear  moduli  of  the  patch 
Transverse  shear  moduli  of  the  sheet 
Shear  modulus  of  adhesive 

Inplane  displacements  of  the  patch  and  sheet  respectively 
Transverse  displacement  of  the  structure 
Axes  of  orthotropy  of  the  patch 

Angle  measured  from  local  axis  system  to  the  material  symmetry  axis  of  the 
patch 

Axes  of  orthotropy  of  the  sheet 

Angle  measured  from  local  axis  system  to  the  material  symmetry  axis  of 
the  sheet 

Element  stiffness  matrix  for  the  adhesive 

Mode  I  stress  intensity  factor  before  and  after  patching,  respectively 
Stress  concentration  factor  before  and  after  patching,  respectively 
Vector  containing  the  degrees  of  freedom  of  a  bonded  element 
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1.  INTRODUCTION 


Advanced  composite  materials  are  finding  increasing  application  in  aerospace  structures. 

In  practice  these  materials  may  contain  defects  such  as  voids,  splits,  delaminations  and 
cracks  produced  by  either  moulding,  machining,  fatigue  or  impact  damage  [1,2].  These  defects 
may  act  as  stress  raisers  which  could  precipitate  failure,  depending  on  the  loading,  geometry 
and  the  “toughness”  of  the  material.  As  a  result  attention  has  recently  been  focused  on  developing 
methods  to  predict  the  residual  strength  of  notched  composites  [3, 4]  and  on  developing  repair 
procedures  [5,6],  In  the  latter  area  the  Aeronautical  Research  Laboratories,  Australia,  have 
pioneered  the  use  of  adhesively  bonded  fibre  reinforced  plastic  patches  to  repair  flaws  in  metallic 
components.  This  procedure  has  been  successfully  used  in  the  repair  of  stress  corrosion  cracks 
in  the  wings  of  Hercules  aircraft,  fatigue  cracks  in  the  landing  wheels  of  Macchi  aircraft  and  in 
the  lower  wing  skins  of  Mirage  111  aircraft  [7,8,9].  The  present  paper  is  concerned  with  the 
extension  of  the  analytical  methods  developed  by  the  authors  for  the  analysis  of  the  repairs  to 
thin  metal  structures  (e.g.  wing  skins)  to  cover  the  case  when  the  structure  is  made  of  a  fibre 
composite  material.  In  this  approach  particular  attention  is  paid  to  the  adhesive  which  bonds 
the  patch  to  the  skin  and  allowance  is  made  for  the  shear  deformation  in  the  skin,  adhesive 
and  patch.  It  is  assumed  that  the  skin  is  made  as  an  orthotropic  laminate  and  that,  in  the  case 
of  a  composite  patch,  the  patch  is  also  an  orthotropic  laminate. 

As  illustrative  examples  the  repairs  of  a  crack  in  a  graphite/epoxy  panel,  and  of  a  hole 
whose  diameter  is  the  same  lenght  as  the  totu '  crack  length  in  the  previous  case  are  studied. 
As  the  repair  material  titanium  sheet,  boron/epoxy  and  graphite/epoxy  laminates  are  considered. 
Various  lay  ups  in  the  damaged  graphite/epoxy  panel  are  investigated.  As  a  result  of  this  analysis 
the  use  of  an  undirectional  laminate  is  recommended  for  repairing  cracks  whilst  the  use  of  titanium 
or  a  quasi-isotropic  lay  up  of  boron/epoxy  is  recommended  for  repairing  holes  in  composites. 


2.  METHOD  OF  ANALYSIS 

Let  us  begin  by  considering  a  thin  composite  patch  which  is  bonded  to  a  thin  sheet  of  fibre 
composite  material.  The  x  and  y  axes  are  taken  in  a  plane  parallel  to  the  midsurface  of  the  sheet 
with  the  z  axis  in  the  thickness  direction.  Under  in-plane  or  transverse  loading,  shear  stresses 
will  be  developed  in  the  adhesive  bond  and  it  is  reasonable  to  assume  that  these  will  be  continuous 
across  the  adhesive-sheet  interface  as  well  as  across  the  adhesive  patch  interface.  Furthermore 
these  shear  stresses,  r«  and  ryi,  are  zero  at  a  free  surface  or  at  a  plane  of  symmetry,  and  it  is 
reasonable  to  assume  that,  since  the  patch  and  the  sheet  are  thin,  these  stresses  vary  linearly 
with  thickness  in  the  patch  and  the  sheet.  With  these  assumptions  the  distribution  of  the  shear 
stresses  r„  and  ryt  is  found  to  be  as  in  [10]:, 


=  f(z)  T,i 

Tv*=/(z)T»r  (0 


where 


/(z)  =  2  z/r,  for  0  ^  z  ^  f,/z 
=  1  for  f«/2  ^  z  <  t,!2  +  «• 

=  (r«/2  +  /»  -f-  to  —  z)//o  for  t»/2  +  f*  <  z  ^  f#/2  +  to  4-  /«  (2) 

and  where  r„,  r,y  are  the  shear  stresses  in  the  adhesive  and  are  assumed  to  be  constant  through 
the  thickness  of  the  adhesive.  The  z  =  0  plane  is  taken  at  the  mid-surface  of  the  sheet  for  a 
doubly  reinforced  sheet  and  at  the  lower  surface  of  the  sheet  for  a  singly  reinforced  sheet.  Here 
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/«,  u  are  the  thicknesses  of  the  patch  and  adhesive  respectively  while  t,  is  the  sheet  thickness  in 
the  doubly  reinforced  case  and  is  twice  the  sheet  thick -&ss  in  the  singly  reinforced  case.  This 
stress  distribution  is  shown  in  Figure  1. 

Let  us  consider  the  case  of  a  doubly  reinforced  sheet.  If  the  x  y  axes  are  at  an  angle  8  to 
the  axes  of  orthotropy,  denoted  by  x'  y\  for  the  composite  patch  and  at  an  angle  g  to  the  axes 
of  orthotropy,  denoted  by  x"  y’,  for  the  composite  sheet  then,  in  the  patch 

r xi  —  C/'i3  (y«  cos2  8  +  yw  sin  8  cos  8)  —  G'aa  (y»*  sin  8  cos  8  —  y«  sin*  8)  (3) 

Tt,  =  G' is  (y«  sin  8  cos  8  -f  yxl  sin*  8)  +  G' 23  (y»*  cos*  8  —  yzt  sin  8  cos  8)  (4) 


where  G' is  and  G' 33  are  the  interlaminar  shear  moduli  of  the  patch  in  the  x'  y  axes  system.  Here 


3w  iu  iw  Sv 

yxt  =  —  +  — ;  yVt  =  —  +  — 

ax  az  ay  az 


(5) 


where  u,  v,  w  are  displacements  in  the  x,  y  and  z  directions  respectively. 

Substituting  for  rx*,and  tx,  as  given  by  equation  (1)  and  for  yxt  and  yw  as  given  by  equation 
(5)  into  equations  (3)  and  (4)  yields,  after  a  slight  rearrangement  of  terms, 


where 


iu  ,  Jv  ,  „  ,  Sw  ,  hw 

r~/x  +  =/(z)T«v  —  /tr - 

Sz  <>z  ax  a>> 


<>« ,  .  av  .  .aw  aw 

r-  /a  +  r~/i  =/(z)  -  hz - /tr~ 

az  az  ax  ay 


(6) 

(7) 


f\  =  (G'13  —  G'23)  sin  8  cos  0 
fz  =  G' 23  cos’  8  +  G'13  sin*  8 

fa  =  G'  13  cos*  8  +  G'33  sin*  8  (8) 

A  similar  system  of  partial  differential  equations  may  be  obtained  for  the  u  and  v  displacements 
in  the  sheet  viz: 


,  a u  av  „  ,  ,  aw  ,  aw 

hi  --  +  >»2  —  —  /(z)  T,y  —  hi- - hz—  (9) 

az  az  ax  ay 

,  a u  av  ,  ,  aw  ,  aw 

ha  —  +  hi  —  =  /(z) t,x  -  ha- - hi  — 

a  z  az  ax  ay 

where 

At  =  (Gia  -  Gz 3)  sin  e  cos  6 
ha  =  C23  cos*  ^  +  613  sin*  g 
ha  =  G 13  cos’  g  +  Gza  sin*  g 

and  where  G13,  G23  are  the  interlaminar  shear  moduli  of  the  sheet  in  the  x'. 

Similarly  in  the  adhesive  layer  we  obtain 


(II) 

y'  axes  system. 


aw  _  aw 

—  =  G«  T„  -  - 

az  ax 


(12) 


av 

az 


=  Ga  t,w  — 


a  w 
ay 


03) 


If  we  now  assume  that  the  vertical  deformation  w  is  independent  of  z,  which  is  consistent  with 
classical  plate  theory  and  Mindlin  plate  theory,  and  that  the  u  and  v  deformations  are  continuous 
across  the  sheet-adhesive  interface  and  the  adhesive-patch  interface  then  equations  (6),  (7),  (9), 
(10),  (12)  and  (13)  may  be  solved  to  give  u  and  v  in  terms  of  w  and  the  shear  stresses  in  the 
adhesive.  Full  details  of  the  solution  process  for  the  case  when  w  =  0  and  the  sheet  is  metallic 
is  given  in  reference  [10].  From  this  solution  it  is  found  that  the  shear  stresses  r,x  and  r,w  may  be 
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DISTRIBUTION  OF  SHEAR  STRESS  THROUGH 
THICKNESS  OF  PATCHED  SHEET 


FIG.  1  AXIS  SYSTEM  IN  PATCH 


related  to  the  displacements  at  the  midsurface  of  the  patch,  which  we  will  denote  by  u0,  vo 
and  w,  and  the  displacements  at  the  midsurface  of  the  sheet,  which  we  will  denote  by  v,  and 
H\  by  the  following  expression: 


Tsx 


Ta* 


=  (uo  -  U, 

=  (uo  -  u. 


+fi  ^zjh/h  + 


+fi 


lx 

l 


^  »’0  -  V, 

/i//a  +  ^v0  —  v. 


+/7  -)m 


,  Iw 

+  fi  r 


0 


hlfi 


(15) 


where  for  convenience  we  have  denoted 


A  =  3/o/i/8  (/e/a  ~/x*)  +  t,  Ai/4  (A2A3  -  A,*) 
h  —  UlGa  +  3tohl%  (hh  —  /1*)  +  U  A2/4  (A2A3  —  A11) 

/•  =  lalGa  +  3/0/3/8  1/2/3  —  /l1)  +  I,  As/4  (A2A3  —  Ai*) 

/?  =  (la  +  1,1 2  +  to/2) 

/s  =  hh  —  A*  (l  6) 

For  a  singly  reinforced  sheet,  i.e.  a  patch  on  one  side  only,  the  term  t,l 4  in  the  expressions 
for  /4,/s  and  /«  is  replaced  by  3r,/8. 

In  the  specific  case  when  bending  effects  are  negligible  and  when  the  sheet  is  isotropic  the 
analysis  reduces  to  that  given  by  the  authors  in  [10].  Furthermore  in  the  case  when  G'is  -  G\ 3 
(  =  Go)  and  G23  =  G 13  (  =  Gt),  which  occurs  when  the  sheet  and  the  patch  are  transversely 
isotropic  we  have 

A  ~A  —  At  =  0,  f2  —  h  —  Go,  h2  —  ha  =  G,,  /s  =/s  =  /o/Ga  +  t,UG,  +  3tolsGo  (17) 
so  that 


T*v 


'  ■  ("°  -  +  ('"  +  I +  f)s)/(s  +  Sf, +  4<0 

-=  ^'0  —  V,  +  +  —  +  2^ G/4G»^ 


(18) 

(19) 


Indeed  in  the  special  case  when  Ju’/Sx  =  aw/a^  =  0  equations  (18)  and  (19)  coincide  with  the 
expression  given,  without  detailed  proof,  in  [12],  Furthermore,  inspecting  equations  (14)  and 
(15)  we  see  that  in  order  for  t,z,  t,v  to  be  zero  we  require 


vo  —  v,  = 


uo  —  u,  = 


(20) 

(21) 


which  is  in  agreement  with  classical  plate  theory. 


3,  ELEMENT  STIFFNESS  MATRIX 


Having  thus  obtained  the  relationship  between  the  adhesive  shear  stress  and  the  displace¬ 
ments  in  the  patch  and  the  sheet  it  is  a  relatively  simple  task  to  formulate  the  stiffness  matrix 
for  an  “adhesive”  element.  We  begin  by  adopting  the  standard  assumption  governing  the  sheet 
and  patch  [10,  12]  viz: 

The  finite  element  model  for  the  sheet  and  the  patch  assumes  a  state  of  plane  stress. 

With  this  assumption  it  is  first  necessary  to  determine  the  strain  energy  due  to  the  shear 
deformation  in  an  element  of  the  structure  viz: 


m 


Txt  y xt  -I-  t„2  Ym  )  dzdxdy 


(22) 
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Here  the  z  intergration  is  over  the  total  thickness  of  the  sheet,  the  adhesive  and  the  patch  while 
the  x  and  _v  integration  is  over  the  area  of  the  element.  Although  in  the  following  analysis  we 
will  only  consider  an  element  which  is  triangular  in  planform  (see  Fig.  2)  the  formulation  may 
easily  be  extended  so  as  to  allow  the  element  to  be  of  any  standard  shape. 


FIG.  2  ELEMENT  OF  ADHESIVE 


Using  the  same  notation  as  in  the  book  by  Zenkiewicz  [13]  we  define  a  vector  f  such  that 

fr  =  (u0,  v0,  us,  0X,  6y) 

=  (Ni)T  (23) 

where  as  explained  in  [13]  page  18,  the  components  of  N  are  generalised  functions  of  position 
and 


5r  -  (6f,  6[,  bTJ 

where 

6,r  =  (Uoi,  *01,  U,t,  1st,  H(,  0zi,  Oyt) 

Here,  as  in  [13],  0X  =  —  iw/iy,  0¥  =  Zw/ix. 

The  strain  vector  y  may  after  some  manipulation  be  expressed  as 

Yr  =  (y„-  =  (°T)r 

where 

*T  =  (r«-  T,y) 

and  the  matrix  D,  which  depends  on  the  thickness  coordinate  r,  is  given  by: 


D  = 


(Af  -  h2hy) 


(■:u) 


(24) 

(25) 

(26) 

(27) 
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in  the  patch 


1  /10V 

'c.^oi/1 


in  the  adhesive 

f(z)  (- 


(28) 


i 

A  h  -hi 


in  the  sheet 


(/?  ~fih)\ 
and  where  f(z)  is  as  defined  in  equation  (I). 

Making  use  of  equations  (14)  and  (15)  and  (23)  we  find  that  x  may  be  written  as 

x  =  ANh 


(29) 


where 

.  _  f  fb'f*'  ~ft  ~~f*'  0’  1 

~  h  U.  fy  -/«.  -/s.  0.  -hhJJi  J 

Substituting  for  x,  as  given  in  (29)  into  equation  (26)  gives 

Y  =  DANS 


(30) 


(31) 


If  we  now  substitute  equations  (1),  (29),  and  (31)  into  the  strain  energy  expression  we  obtain 


V  = 


f(z)  DdzJ  ANh  dxdy 


(32) 


The  stiffness  matrix  Kr  may  now  be  obtained,  as  outlined  in  (13),  by  differentiating  V  with 
respect  to  S  and  is  given  by 


r  =  ||  (/<N)7||  f(z)  Ddz^j  AN  dxdy 


(33) 


As  before  the  double  integration  is  over  the  area  of  the  element  while  the  integration  with  respect 
to  z  is  over  the  thickness  of  the  structure. 

This  formulation  looks  quite  complex  but  in  the  case  when  the  interlaminar  shear  moduli 
of  the  patch  (and  the  sheet)  are  equal  it  simplifies  considerably.  For  example,  let  the  structure  be 
patched  on  both  sides  and  let  G'n  =  G'23  (=  G0)  and  Cu  =  G23  (=  G,).  This  gives 

K'  =  2  [(tJGa  +  r0/3C0  +  /,/6G^||  (A N)T  ^  AN  dxdy  (34) 

where  A  has  simplified  to 


A  = 


1,0,  —I,  0, 0, 0,/7  [a  t  3r0\ 

.0,  1,0,  -1,0,  -/„  0//  \Ga  4 Gs  8 G0) 


and  as  before  /7  =  ta  +  tj 2  +  IJ2.  It  is  important  to  note  that  as  the  adhesive  thickness  tends 
to  zero  (i.e.  ta  -*  0)  the  stiffness  matrix  Ke  tends  to  a  constant  and  non  zero  value.  This  should 
enable  the  present  approach  to  be  used  to  model  internal  delaminations  as  well  as  the  repair 
problems  which  will  be  considered  here. 

This  approach  coincides  with  the  analysis  presented  in  [10]  when  iw/ix  =  Zw/Zy  =  0  and 
the  sheet  is  isotropic.  Furthermore  when  Zw/ix  =  Zw/Zy  =  0  but  the  sheet  is  a  composite 
laminate  the  functional  form  of  the  shear  stresses  tsx  and  r„,  as  given  by  equations  (14)  and 
( i 5),  is  similar  to  that  given  in  [10].  The  difference  is  that  the  expressions  for /4,/5,  and /6  given 
in  the  present  paper  are  more  complex  than  those  given  in  [10]  for  an  isotropic  sheet. 


4.  ILLUSTRATIVE  EXAMPLES 

Let  us  now  consider  the  repair  of  damaged  fibre  composite  panels.  The  basic  ply  is  a  graphite/ 
epoxy  (Narmco  Thomel  T300/5208  with  the  following  material  properties  £,  —  141  0  x  I03 
MPa ,  ul2  =  0-31,  G,j  =  518  x  10 *MPa,  £|/£2  =  14-96.  Since  it  is  often  stated  that  com¬ 
posites  are  relatively  notch  insensitive  two  particular  kinds  of  damage  were  considered.  The 
first  is  a  centrally  located  crack  38  - 1  mm  long  in  a  graphite  epoxy  panel  of  dimensions  508  mm  x 
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635  mm  2-29  mm.  The  second  kind  of  damage  is  a  centrally  located  hole  with  a  diameter 
of  38-1  mm  in  the  same  panel.  Various  panel  lay  ups  were  considered  but  to  enable  a  simple 
comparison  of  the  numerical  results  the  panel  thickness  was  taken  as  2-29  mm  for  each  lay  up. 
The  moduli  of  the  panel  were  taken  to  be  as  for  a  symmetric  lay  up  of  either 

(1)  (0  ;  45.90), 

(2)  (0  90), 

(3)  (0, ...  45), 

(4)  (x  45), 

(5)  (0), 

These  laminates  were  chosen  because  of  their  similarity  to  skin  lay  ups  in  existing  or  planned 
military  aircraft.  In  each  case  the  panels  were  assumed  to  be  subjected  to  a  uniform  tensile  stress 
at  these  edges.  Because  of  the  symmetrical  nature  of  these  problems  only  one  half  of  the  structure 
was  analysed. 

The  finite  element  mesh  for  the  cracked  panel  consisted  of  147  constant  strain  triangles, 
142  quadrilateral  elements  and  a  special  crack  element,  (see  reference  [10].) 

The  mesh  for  the  panel  containing  a  hole  consisted  of  approximately  the  same  number  of 
elements. 

To  these  models  was  added  a  finite  element  model  of  the  various  bonded  repairs  (i.e.  patches). 
A  patch  was  assumed  to  be  placed  on  both  sides  of  the  pane).  Each  patch  is  stepped  in  thickness 
and  has  dimensions  152  mm  •  50-8  mm  with  a  maximum  thickness  of  0  762  mm.  which  occurs 
in  the  middle  of  the  patch;  see  Figure  3.  The  patches  considered  were; 

(a)  Titanium 

(b)  A  unidirectional  graphite  epoxy  patch  with  the  fibres  in  the  direction  of  the  load. 

(r)  A  graphite/'epoxy  laminate  with  the  same  moduli  as  the  sheet 

( d )  A  quasi-isotropic  (e  g.  0  ;  45/90)  boron  epoxy  laminate. 

The  finite  element  model  of  these  patches  consisted  of  150  constant  strain  triangles  and  150  of 
the  adhesive  elements  described  above.  In  this  problem  the  analysis  simplified  considerably 
since  the  sheet  was  patched  on  both  sides  and  the  load  was  acting  in  the  plane  of  the  midsurfacp 
of  the  sheet.  As  a  result  the  structure  did  not  experience  bending  effects,  i.e. 


iw 

w 

S.x 


(36) 


This  in  turn  simplified  the  Formulation  of  the  adhesive  stiffness  matrix,  since  it  was  now  only- 
necessary  to  specify  the  way  in  which  the  wand  r  displacements  varied  in  the  adhesive  element. 
In  this  paper  we  follow  the  approach  given  in  [10]  using  a  triangular  element  and  taking 

u  ((a/  t-  bt  x  f  G.r)w(  ;  («;  t  bj\  ■  cjy)  uj  (a™  h,„  x  ■  Cm.v)  um)  2A  (37) 

and  a  similar  expression  for  v  where  as  in  [13] 

at  X)Ym  —  XmVj.  bt  *  Y,„  -  Yj.  Ci  Xm  Xj  (38) 

and  flj,  <7m,  bj.  bm.  cj  and  r,„  are  obtained  by  a  cyclic  permutation  of  /.  /.  m.  Here  (a,,  v,  ).  (xj.Yj) 
and  (xm-  Xm)  are  the  coordinates  of  the  corners  of  the  element,  A  is  the  area  and  wi.  v,  are  the 
displacements  at  the  ith  node  (see  Fig.  2). 

The  effect  that  patches  have  on  the  damaged  panels  is  described  below  where  we  first  consider 
the  repair  to  the  cracked  panel  using  either  patch  (a)  or  patch  (6).  In  Table  I  we  see  the  ratio  of 
the  stress  intensity  factors  KtpIKiu.  where  Ktp  and  Ktu  are  the  values  of  the  stress  intensity  factor 
after  and  before  patching  respectively.  The  values  obtained  for  Ktu  coincided  with  those  pre¬ 
viously  published  and  which  differ  from  the  value  obtained  for  an  isotropic  sheet  by  only  a  few 
percent. 
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FIG.  3  PATCHED  PANEL 


Table  1.  Ratios  of  stress  intensity  factors 

A] p  Kiu 


Sheet  Material 

Repair  Material 

Unidir.  graphite  epoxy 

Titanium 

(0  :  45/90), 

016 

0  14 

(090), 

0- 14 

0  14 

<0  :  45), 

0-17 

016 

(.1  45). 

0-11 

Oil 

(0), 

0-18 

0-17 

From  this  table  we  see  that  as  a  repair  to  cracked  graphite  epoxy  sheets  the  unidirectional 
graphite  laminate,  with  the  fibres  perpendicular  to  the  crack,  viz  in  the  direction  of  the  load, 
and  the  titanium  are  almost  equally  effective  in  reducing  the  stress  intensity  factor  at  the  crack 
tip.  This  is  consistent  with  the  results  given  in  [II]  which  deals  with  the  repair  of  a  crack  of  the 
same  length  in  an  aluminium  panel  of  the  same  dimensions  as  considered  here.  Indeed  in  [II] 
the  authors  show  that  a  slight  increase  in  the  stiffness  of  the  patch  does  not  greatly  reduce  the 
stress  intensity  factor.  Nevertheless  at  first  glance  it  is  surprising  that  the  titanium,  which  has 
a  Young’s  modulus  of  106-7  ■  103  MPa.  is  as  effective  as  the  unidirectional  graphite  epoxy 
laminate  with  Fu  141-0  •  103  MPa.  The  explanation  of  this  phenomenon  is  partly  due  to 
the  fact  that,  as  indicated  above,  the  rate  of  reduction  of  the  stress  intensity  factor  with  increasing 
stiffness  is  small  and  partly  due  to  the  fact  that  the  transverse  shear  modulus  for  the  titanium 
(Gig  =  Ci 3  ~  40-5  103  MPa)  is  much  greater  than  that  of  the  graphite  epoxy  (G< 3  =  C]3 

=  5-2  x  103MPa).  In  titanium  the  crack  opening  displacement  is  resisted  by  virtually  the  entire 
thickness  of  the  titanium,  whereas  in  graphite  the  low  value  of  its  transverse  shear  modulus 
results  in  the  plies  being  unequally  stressed  with  the  plies  closest  to  the  surface  providing  the 
most  resistance  to  the  opening  of  the  crack. 

Although  here  we  have  confined  our  attention  to  unidirectional  graphite  patches  there  are 
many  reasons  for  using  a  uni-directional  boron  patch  to  repair  cracks.  One  of  the  main  reasons 
is  that  standard  eddy  current  inspection  procedures  may  be  used  to  inspect  for  crack  growth 
under  a  boron  patch.  For  further  details  concerning  the  advantages  of  boron  epoxy  over  carbon 
epoxy  laminates  as  repair  material  see  [7], 

Although  we  have  seen  that  both  repairs  give  rise  to  approximately  the  same  reduction  in 
the  stress  intensity  factors  there  are  several  reasons  why  a  unidirectional  repair  should  be  used 
in  preference  to  a  titanium  patch.  Indeed  the  reasons  for  recommending  use  of  a  uni-directional 
laminate  to  repair  cracked  composite  sheets  are  the  same  as  for  recommending  the  use  of  a 
unidirectional  laminate  to  repair  cracked  metallic  sheets  and  are  given  in  [7.  8).  These  are  best 
illustrated  by  considering  the  recent  repair  to  cracks  in  the  lower  wing  skin  of  Mirage  III  aircraft 
in  service  with  the  Royal  Australian  Air  Force;  see  reference  [8],  This  repair  was  a  unidirectional 
boron/epoxy  laminate  with  the  fibres  perpendicular  to  the  crack  and  was  designed  by  the  authors. 
In  this  case  the  crack  lay  in  vicinity  of  the  spar  and  the  root  rib  intersection.  The  use  of  a  titanium 
patch,  with  its  parasitic  stiffness  in  the  direction  parallel  to  the  crack,  could  easily  have  changed 
the  strain  distribution  in  the  spar,  which  was  itself  stress  critical.  However  by  using  the  unidirec¬ 
tional  boron  laminate  the  strain  in  the  spar,  with  a  cracked  but  patched  wing  skin,  was  restored 
to  the  value  of  the  strain  in  a  spar  with  an  uncracked  wing  skin;  see  reference  [8], 

Let  us  now  turn  our  attention  to  the  repair  of  holes  in  graphile/epoxy  laminates.  In  Tables  2 
and  3  are  given  the  ratio  of  the  stress  concentration  factors  AYr/AVr  and  the  ratio  of  the  maximum 
fibre  stress  to  the  applied  stress  <?plo.  Here  Ktp  and  Krr  are  the  stress  concentration  factors  in 
the  sheet  after  and  before  patching  respectively.  The  values  obtained  for  Krt<  coincided  with  those 
previously  given  in  [3]. 


Table  2.  Ratio  of  Stress  Concentration  Factors  KtpIKti 


Panel  laminate 


Patch 

0/d  45 

r  45 

0/90.  t  45 

0/90 

0 

a 

0-32 

0-32 

0-33 

0-32 

0-26 

b 

0-30 

0-28 

0-32 

0-26 

0-26 

c 

0-43 

0-54 

0-46 

0  40 

0-26 

d 

0-36 

0-32 

0-32 

0-37 

0  31 

Table  3.  Ratio  of  Patch  Stress  to  Sheet  stress  ap  'a 


Panel  laminate 
Patch 

0/.»  45 

;  45 

0  90/  ;  45 

0  90 

0 

a 

3  05 

4-2 

2  69 

2-25 

2-32 

b 

3-94 

5-4 

3-46 

2-86 

3  08 

c 

2-43 

1  -73 

1  91 

2  02 

3  08 

d 

2-73 

3-82 

3-46 

201 

212 

From  Tables  2  and  3  we  see  that  when  repairing  holes  in  composite  sheets  under  uniaxial 
loading  it  is  best  to  use  a  unidirectional  laminate  with  the  fibres  in  the  direction  of  the  load 
Indeed  one  may  generalize  this  and  say  that  if,  under  an  arbitrary  system  of  loads,  one  knows 
the  direction  of  the  principal  stress  which  is  primarily  responsible  for  failure  then  it  is  best  to 
use  as  a  repair  a  unidirectional  laminate  with  the  fibres  in  the  direction  of  the  principal  stress. 
In  general  however  one  may  not  know  this  direction  and  in  this  case  the  use  of  a  unidirectional 
laminate  as  a  repair  could  be  unconservative.  Indeed  in  this  situation  the  use  of  a  bonded 
titanium,  or  a  quasi-isotropic  boron  patch  as  a  repair  is  quite  common  [5],  Furthermore  consulting 
Table  2  we  see  that  the  titanium  and  the  quasi-isotropic  boron  patches  give  approximately  the 
same  ratio  of  the  stress  concentration  factors  which  is  significantly  lower  than  that  obtained 
using  a  patch  with  the  same  lay  up  as  the  skin  material. 

As  a  result  since  bonded  repairs  often  take  place  either  in  the  field  [5,  8]  where  the  direction 
of  the  principal  stress,  may  not  be  available,  or  in  a  repair  depot  the  recommended  repairs  to 
holes  in  composite  sheets  is  either  a  bonded  titanium  or  a  quasi-isotropic  boron  patch.  On  the 
other  hand  if  the  panel  contains  a  distinct  crack  the  recommended  repair  is  a  unidirectional 
laminate  with  the  fibres  perpendicular  to  the  crack. 


5.  CONCLUSION 

An  advanced  finite  element  method  for  analysing  the  repair  of  damaged  composite  laminates 
has  been  developed.  Subsequent  analysis  has  shown  that  either  a  titanium  or  a  quasi-isotropic 
boron  repair  is  equally  suited  to  the  repair  of  holes  while  a  unidirectional  laminate,  with  its 
fibres  perpendicular  to  the  crack,  is  best  suited  to  the  repair  of  cracks  in  composite  laminates. 
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